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Abstract 

This paper addresses the problem of bearing-based network localization, which aims to localize all the nodes in a static network 
given the locations of a subset of nodes termed anchors and inter-node bearings measured in a common reference frame. The 
contributions of the paper are twofold. Firstly, we propose necessary and sufficient conditions for network localizability with 
both algebraic and rigidity theoretic interpretations. The analysis of the localizability heavily relies on the recently developed 
bearing rigidity theory and a special matrix termed the bearing Laplacian. Secondly, we propose a linear distributed protocol 
for bearing-based network localization. The protocol can globally localize a network if and only if the network is localizable. The 
sensitivity of the protocol to constant measurement errors is also analyzed. One novelty of this work is that the localizability 
analysis and localization protocol are applicable to networks in arbitrary dimensional spaces. 
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1 Introduction 

Distributed localization of sensor networks is a core problem in many multi-agent coordination tasks. Network 
localizability and distributed protocols are two fundamental problems for any network localization problems. Network 
localizability characterizes whether or not a network can be possibly localized given the anchor locations and inter¬ 
neighbor relative measurements, whereas distributed protocols are used for localizing the network in a distributed 
manner if the network is localizable. According to the types of the relative measurements used for localization, the 
existing works can be divided into three classes: distance-based, bearing-based, and position-based. Distance-based 
network localization has been studied extensively so far (see [1-4] and the references therein). The analysis of the 
localizability in distance-based network localization relies heavily on the distance rigidity theory. It has been shown 
that a network in an n-dimensional space can be uniquely localized if the network is globally rigid and has at 
least n + \ anchors in a general position [1]. More recently, bearing-based network localization has also attracted 
extensive research attention [5-11]. The analysis of the localizability in bearing-based network localization relies on 
the analogous bearing rigidity theory [12-15]. Finally, position-based network localization, where the inter-neighbor 
distance and local bearing measurements are used together for network localization, has been studied in [16] by using 
a complex graph Laplacian. 

Although bearing-based network localization has been studied by many researchers, the two fundamental problems, 
network localizability and distributed protocols, have not yet been fully explored. It was shown in [7-10] that a 
network is localizable when the network is bearing rigid and has at least two anchors. This condition is, however, 
sufficient but not necessary when the number of anchors is greater than two [10, Cor 10]. A necessary and sufficient 
condition for network localizability was proposed in [10, Thm 15] based on the notion of a stiffness matrix. This 
condition is, however, applicable only to networks in two-dimensional spaces. In fact, the localizability of a network 
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is jointly determined by many factors such as its topological and Euclidean structure, as well as the selection of the 
anchors. The relationship between the localizability and these factors have not been fully understood yet up to now. 
Moreover, the existing bearing-based localization protocols are mainly applicable to networks in two-dimensional 
ambient spaces [7-10]. General results of localizability or distributed protocols for bearing-based network localization 
in three and higher dimensional spaces are still lacking. 

This paper studies the localizability and distributed protocols for bearing-based network localization in arbitrary 
dimensional spaces. It is assumed that the anchors’ locations and inter-neighbor bearings measured in a global 
reference frame are already given. The main contributions of this work are summarized below. 

(a) We first show that the bearing-based network localization problem can be formulated as a linear least-squares 
optimization problem. A special matrix termed the bearing Laplacian, which can be viewed as a matrix-weighted 
graph Laplacian, emerges as a key part in the least-squares formulation and plays important roles in the subse¬ 
quent analysis. 

(b) Based on the least-squares formulation, we propose necessary and sufficient conditions for network localizability 
with both algebraic and rigidity theoretic interpretations. These conditions not only provide numerical ways to 
examine the localizability of a given network but also provide intuitions on what a localizable network looks like. 

(c) We then propose a distributed linear localization protocol. It is proved that the protocol can globally localize a 
network if and only if the network is localizable. The sensitivity of the protocol to constant measurement errors 
is also analyzed. 

Finally, it is worth noting that the localizability analysis presented in this paper is independent to whether the 
sensing graph is directed or undirected. The convergence analysis of the proposed localization protocol, however, 
relies on the assumption of undirected sensing graphs. The convergence analysis of the protocol in the directed case 
is considered in [17]. 

The rest of the paper is organized as follows. Section 2 presents the linear least-squares formulation of the bearing- 
based network localization problem. Section 3 analyzes the properties of the bearing Laplacian and its connection to 
the bearing rigidity theory. Section 4 presents necessary and sufficient conditions for network localizability. Section 5 
proposes and analyzes a linear distributed localization protocol. Conclusions are drawn in Section 6. 


Notations: Given Ai € for i = 1, ... ,n, denote diag(Ai) = blkdiag{Ai,..., A„} € Let jj • jj be the 

Euclidian norm of a vector or the spectral norm of a matrix, and 0 be the Kronecker product. Denote Id € 
as the identity matrix, and 1^ = [1, ■. ■, I]"*" G K'’*. Let Null(-) and Range(-) be the null space and range space of a 
matrix, respectively. 


2 Problem Formulation of Bearing-Based Network Localization 

In this section, the problem of bearing-based network localization is formally stated and then formulated as a linear 
least-squares problem. Central to this problem is the notion of localizability, which is formally defined here. 

2.1 Problem Statement 


Consider a network of n stationary nodes in {n > 2 and d> 2). Assume no two nodes are collocated. Let pi G M'’* 
be the location of node i (i = 1,... ,n). Define the edge vector and the bearing between nodes i and j as 


^ij — Pj Pi 1 



The unit vector gij represents the relative bearing of pj with respect to pi. Note = —eji and gij = —gji. 
Suppose the locations of Ua anchor nodes are already given and the locations of the remaining n/ follower nodes 
are to be estimated {ua + Uf = n). Denote Va = {1, ■ • ■ ,na}, V/ = {ua -I- 1,... ,n}, and V = Va U V/. Denote 


Pa=[pJ,...,plJ^ G 


\Pf = [p^ 


ria + l ’ ■ 


,pI? e 


rp T 


, and p = [pj,p]] 


Suppose each node has the bearing-only sensing capabilities. The sensing topology of the network defines a graph 
Q — iy,S) where 5 C V x V. Denote (i,j) as the directed edge with node i as the tail and node j as the head. The 
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Fig. 1. An illustration of the notion of localizability. Black dots represent the anchors and white dots for the followers. Suppose 
the true network is (a). The networks in (a) and (b) both satisfy the nonlinear equations in (1). The networks in (a), (b), and 
(c) all satisfy the linear equations in (2). 


directed edge (z, j) G £ indicates that node i can “see” node j; that is node i can measure the relative bearings 
of node j. Node j is called the neighbor of node i if (z, j) G £, and A/i = {j G V\{i,j) G £} is the neighbourhood 
of node i. We assume a global orientation that can be sensed by all the nodes, and thus all measured bearings 
can be expressed with respect to this common orientation. The global orientation means a common north for the 
two-dimensional space, and a common north-east-down reference for the three-dimensional space. Finally, let Q{p) 
denote the network that is the graph Q with each vertex z G V mapped to the point pi. 


The problem of bearing-based network localization is formally stated below. 


Problem 1 (Bearing-Based Network Localization) Consider a network G{p) in the bearing-based network 
localization problem is to determine the locations of the follower nodes, {pijigV/; given the inter-neighbor bearings, 
{gij}(i,j)^e! nnd the locations of the anchor nodes, {pi}igVa- Mathematically, the problem is to retrieve the true 
network location p by solving the system of nonlinear equations. 


Pj - Pt 
\\Pj-P^\\ 
Pz =P^, 


9ij J J ) G 
Vf G Va, 


( 1 ) 


where Pi is the estimated location of node i. 

The true network location is always a solution to the nonlinear equations in (1), but the nonlinear equations may 
admit many other solutions that do not correspond to the true network location. Thus we need to study when the 
true network location is the unique solution to (1), which motivates the following notion. 

Definition 1 (Bearing-Based Network Localizability) A network Q{p) is called bearing-based localizable if 
the true network location p is the unique solution to (1). 

Localizability is a fundamental property of bearing-based networks. A network must be localizable in order to be 
localized with either distributed or centralized protocols. The notion of localizability is illustrated by an example in 
Figure 1. In this example, the network in Figure 1(a) is the true network. The network in Figure 1(b) has the same 
bearings and anchor locations as the true network. As a result, both of the networks in Figure l(a)-(b) are solutions 
to (1) and hence the networks are not localizable by Definition 1. 

For the sake of simplicity, we assume that the graph Q is undirected, which means (z,j) G £ (j,z) G £. If the 

graph is directed, suppose (z, j) G £ but (j,z) ^ £. We can always add the edge (j, z) into £ to convert the directed 
graph to an undirected one. The directed edges {i,j) and (j,z) imply two equations {pj —Pi)/\\Pj — Pi\\ = gij and 
{Pi ~ Pj)l\\Pi ~ Till = 5ii) respectively. The two equations are equivalent because gji = —gij- As a result, adding the 
edge (j,z) does not affect the solutions to (1). 

2.2 Reformulation as a Least-Squares Problem 

In order solve the nonlinear equations in (1), we derive a companion system of linear equations. In this direction, we 
first introduce a useful orthogonal projection operator. For any nonzero vector a; G {d> 2), define the orthogonal 


3 









projection operator P : —>• as 


P{x) ^ Id 


T 

X X 

Ikll Ikll ■ 


For notational simplicity, denote Px = P{x). The matrix Px geometrically projects any vector onto the orthogonal 
compliment of x. It can be easily verified that Pj = Px, P^ = Px, Null(Pa;) = spanja:}, and the eigenvalues of Px 
are {0, 

Consider now the projection matrix, Pg.^ = Id — dijdjj, associated with the bearing gij. By multiplying Pg.^ on both 
sides of the first equation in (1), the nonlinear algebraic problem (1) is converted to a system of linear equations. 


I Pg,, (Pi - Pi) = j) e £, ^2) 

\pi=Pi, Vi e Va- 

System (2) is not equivalent to system (1) in general. But we have the exact relation between (1) and (2) as described 
in the following lemma. 

Lemma 1 Let Xi and X 2 denote the set of all solutions satisfying (1) and (2), respectively. Then 

(a) {p} C Ti C Ta; 

(b) {p} = Xi if and only if {p} = X 2 . 


Proof, (a) Since the true network location p is always a solution to (1) and (2), we know Xi and X2 are nonempty 
and {p} C Xi and {p} C X 2 . Since (2) is obtained by multiplying (1) by Pg.., we know any solution to (1) is also a 
solution to (2), showing Xi C T 2 . 

(b) (Sufficiency) Suppose {p} = X 2 . It then follows from {p} C dfi C ^”2 that {p} = Xi. (Necessity) Suppose 
{p} = Xi . We next prove {p} = X2 by contradiction. Assume p' G X2 and p' 7^ p. Let 6p = p' — p and define 

p" = p + kSp, fc € M. (3) 

We next show that p” G Xi when |fc| is sufficiently small, leading to a contradiction. Since p,p' G X 2 , we know 
p” G X 2 for all /c G M by (3). As a result, for any fc G M and (i,j) G 6, we have Pg..(p" — p'f) = 0 which implies 
either (p" - p")/lb" - p”\\ = or (p” - p'l)/\\p'- - p"|| = -gij. Since p" - p" = (pj - pf) + k(6pj - 6pf) according 
to (3), it is obvious that when |A:| is sufficiently small, the entries of p" — p" have the same signs as those of pj — pi, 
and consequently (p" — p'DfWp'j — p'fW = (pj — Pi)/\\Pj — Pi\\ = 9ij- Note that when any entry of pj — pi is zero, 
the corresponding entry of dpj — 6pi is also zero because Spi — dpj is parallel to pj — Pi- To conclude, p" is another 
solution other than p satisfying (I), which is a contradiction. □ 

Remark 1 The proof of Lemma 1(b) can be illustrated by Figure 1, where the networks (a), (b), and (c) correspond 
to p, p”, and p' in the proof, respectively. 


Lemma I indicates that the true network location p is the unique solution to (1) if and only if p is the unique solution 
to (2). Thus we can study the localizability by analyzing the linear system (2). The linear system of equations in (2) 
can be rewritten as the following linear least-squares problem. 


minimize 
subject to 


J{P) = 11^90b* 

iev jeW 


Pi=Pi, z G V, 


(4) 


Since any minimizer with the objective function as zero is the solution to (2), we now successfully formulate the 
localizability problem as the above least-squares problem. The rest of the paper is dedicated to studying two prop- 
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erties of the least-squares problem. The first is to determine when the true location p is the unique global minimizer 
of (4) (i.e., the network is localizable), and the second is how to obtain p in a distributed manner. 


3 The Bearing Laplacian Matrix 

In this section, we show that a new important matrix, termed bearing Laplacian, emerges in the least-squares 
formulation. The useful properties of the bearing Laplacian that will be used throughout the paper are explored. 

Since the underlying graph Q is undirected, the objective function in (4) can be expressed in a quadratic form, 

J(p) =f'B{g{p))p, 

where B{Q{p)) G and its ijth subblock matrix is 

{ Odxdj I ^ J; (hj) ^ 

-Pgij, (Li)e£, 

J2ke^r.Pg^k, i=j,iGV. 

For notational simplicity, we write B{g{p)) as B in the sequel. The matrix B has a structure reminiscent of the 
weighted graph Laplacian matrix. Since B indicates not only the topology of the network but also the inter-neighbor 
bearings, it is referred to as hearing Laplacian in this paper. 

The bearing Laplacian has an intimate connection to the bearing rigidity properties of the network. Preliminaries 
to the bearing rigidity theory, originally proposed in [15], are given in Appendix A. Here we would like to highlight 
two important notions from this theory. The first is the notion of infinitesimal bearing motions. Loosely speaking, 
infinitesimal bearing motions are motions of the nodes that preserve inter-neighbor bearings. For example, for the 
network in Figure 1(a), the bearings can be preserved when the nodes 3 and 4 move in the horizontal direction 
to the right. A network always has two kinds of trivial infinitesimal bearing motions - they are the translational 
and scaling motions of the entire network. A network is infinitesimally bearing rigid if all its infinitesimal bearing 
motions are trivial. One important property of an infinitesimally bearing rigid network is that its shape can be 
uniquely determined by the inter-neighbor bearings. 

We next give the basic properties of the bearing Laplacian matrix. We also show that the bearing Laplacian matrix 
is a powerful tool for characterizing the bearing rigidity of a network. 

Lemma 2 For a network g{p) with undirected graph Q , the bearing Laplacian B satisfies the following: 

(a) B is symmetric positive semi-definite; 

(b) Rank(;B) < dn — d — 1 and Null(;B) A span{l (g) Ld,p}; 

(c) Rank(,B) = dn — d — 1 and Null(S) = spanjl (g Id,p} if and only if Q{p) is infinitesimally bearing rigid. 


Proof. Assign an arbitrary orientation to each undirected edge and label the edge vectors and bearings for the 
directed edges as {ek}'k=i and {gk}'k=n respectively. Then the bearing Laplacian B can be expressed as B = 
H'^d\ag{Pg,.)H where H = H®Id and H is the incidence matrix of the graplrPHlt further follows from Pg,. — Pj^Pg^ 
that 

B = iL'^diag(Pjj diag(Pg JiL = 
nr n 


^ The incidence matrix H £ Qf an oriented graph is the {0, ±l}-matrix with [H]ki = 1 if vertex i is the head of edge 

k, = —1 if it is the tail, and 0 otherwise. 


5 





Note TZ = diag (||efe|j/d) Rb where Rb is the bearing rigidity matrix (see Lemma 7 in Appendix A). As a result, the 
matrix TZ, and hence B, have exactly the same rank and null space as Rb- Then the results in (b) and (c) follows 
immediately from Lemma 7 and Theorem 7 as given in Appendix A. □ 


Since the nodes in the network are partitioned into anchors and followers, it will be useful to partition the corre¬ 
sponding bearing Laplacian as 


^ _ Baa Baf 
~ Bff_ ’ 

where Baa G Baf = Bj^ G g j^dufXduf^ 

Lemma 3 For any network Q(p) with undirected graph Q, the subblock matrix Bf f is symmetric positive semi-definite 
and satisfies BffPf + BfaPa = 0. 

Proof. For any nonzero x G denote x = G Since ,8 > 0, we have x'^BffX = x^Bx > 0. As 

a result Bjf is positive semi-definite. Since p G Null(8) as suggested by Lemma 2, we have Bp = Q which further 
implies BfaPa + BffPf = 0. □ 


4 Analysis of Network Localizability 

In this section, we analyze the localizability of networks in arbitrary dimensions. We hrst prove two necessary and 
sufficient conditions for network localizability from algebraic and rigidity perspectives, respectively. We then present 
more necessary and/or sufficient conditions which can give more intuition on what localizable networks look like. 

First of all, we derive the optimality condition for the least-squares problem (4). 

Lemma 4 For the least-squares problem (4), any minimizer p*j is also a global minimizer and satisfies 


BffP*f + BfaPa = 0. 


Proof. By substituting pa = Pa into the objective function J{p) = p^Bp, the constrained optimization problem (4) 
can be converted to the unconstrained problem 

min J{pf)= p] BffPf + 2plBafPf + P^BaaPa- (5) 

p f 


Any minimizer must satisfy 'Z/pfJ{pf) = BffPf -\- BfaPa = 0. Now suppose is a minimizer and satisfies Bffp*j -\- 
BfaPa = 0. By comparing with BffPf -\~BfaPa = 0 as shown in Lemma 3, we know p*j = pf -\-x where x G Null(8//). 
Let p* = [Pa: (P/s-nd x = G Since p*f = pf -\- x and Bp = 0, we have J{p*) = {p*)'^Bp* = 

{p -\- xp-B{p -I- a;) = x^Bx = x^BffX = 0. As a result, the objective function equals zero at every minimizer. □ 


The linear equations in (2) hold if and only if the objective function in the least-squares problem (4) is minimized to 
zero; this is a direct consequence of the hrst-order optimality conditions associated with (4). Thus the equivalence 
between (2) and (4) is formally established. We are now ready to present the necessary and sufficient condition for 
localizability. 

Theorem 1 (Algebraic Condition for Localizability) A network Q{p) is localizable if and only if the matrix 
Bff is nonsingular. When the network is localizable, the true locations of the followers can be calculated by pf = 
- Bjj BfaPa- 
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(a) 




(d) 


(e) 


(f) 


Fig. 2. Examples of non-localizable networks. The solid and hollow dots represent the anchors and followers, respectively. The 
networks are not localizable because they have infinitesimal bearing motions that only correspond to the followers (see, for 
example, the red arrows). The networks in (e) and (f) are three-dimensional, and the rest are two-dimensional. 



Fig. 3. Examples of localizable networks. The solid and hollow dots represent the anchors and followers, respectively. The 
networks in (e), (f), and (h) are three-dimensional, and the rest are two-dimensional. 

Proof. By Lemma 4, a network is localizable if and only if the true network location p is the unique minimizer of 
the least-squares problem (4). Since any minimizer must satisfy Bffp*j + BfaPa = 0, it is obvious that the minimizer 
is unique if and only if S// is nonsingular. When Bjf is nonsingular, we have = —BjjBfaPa, whose value equals 
the true location pf according to Lemma 3. □ 


Theorem 1 establishes the equivalence between the localizability and the nonsingularity of Bff. A question that 
immediately follows Theorem 1 is what kind of networks have nonsingular Bff. We next propose a necessary and 
sufficient condition from the bearing rigidity point of view. This rigidity condition is mathematically equivalent to 
the algebraic condition, but it gives more intuition on what localizable networks look like. 

Theorem 2 (Rigidity Condition for Localizability) A network Q (p) is localizable if and only if every infinites¬ 
imal hearing motion involves at least one anchor; that is, for any nonzero infinitesimal bearing motion 


5p = 


dpa 

6pf 


G Null(S), 


the vector Spa corresponding to the anchors must he nonzero. 


Proof. We only need to show that Bff is singular if and only if there exists nonzero Sp G Null(,B) with 6pa = 0. 
{Necessity) Suppose Bff is singular. Then there exists nonzero x G such that BffX = 0. Let 5p = [0, G 
Then 5p^B5p = x^BffX = 0. Hence 5p G Null(H) and 5pa = 0. (Sufficiency) Suppose there exists Sp G Null(H) 
satisfying Spa = 0 and Spf 7 ^ 0. Then SpJBffSpf = Sp'^BSp = 0, which implies that Bff is singular. □ 


The intuition behind Theorem 2 is as follows. Any infinitesimal bearing motion (i.e., bearing-preserved motion) 
would imply multiple false networks that have exactly the same bearings as the true network. Only if the infinitesimal 
bearing motion involves at least one anchor, the false networks can be ruled out as solutions to ( 1 ) since they do not 
satisfy the anchor constraints; otherwise, the false networks cannot be distinguished from the true network. 

Examples are given in Figure 2 and Figure 3 to illustrate Theorem 2. Figure 2 shows examples of non-localizable 
networks. These networks are not localizable because each of them has infinitesimal bearing motions that only involve 
the followers (see those marked by red arrows). Figure 3 shows examples of localizable networks. The networks in 
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Figure 3(a)-(f) are obtained by modifying the networks in Figure 2, which suggests that a non-localizable network 
can be made localizable by adding extra edges or selecting different anchors. It is worth noting that the networks in 
Figure 3(c)-(g) are not infinitesimally bearing rigid yet they are localizable. As a result, infinitesimal bearing rigidity 
is not necessary to guarantee localizability. 


Up to this point, we have presented two necessary and sufficient localizability conditions. One is the algebraic 
condition in Theorem 1 and the other is the rigidity condition in Theorem 2. We next utilize the two conditions to 
examine some specific problems more closely. The first is how many anchors are required to ensure the localizability 
of a network. 

Corollary 1 If a network Q (p) is localizable, then 


dim (Null(, 8 )) 

na > J 


> 1 . 


Proof. Let k = dim(Null(, 8 )) and N G '^dnxk basis matrix of Null(, 8 ) which means Range(A^) = Null(,B). 
Then any nonzero Sp G Null(,B) can be expressed as Sp = Nx, where x G M^,a; ^ 0. Partition N and express 


Nx as Sp = Nx = 


NaX 

Nfx 


where Na G According to Theorem 2, the network is localizable if and 


only if NaX 7 ^ 0,Va; G yf 0. As a result, the matrix Na must have full column rank, which requires Na 

to be a tall matrix with dua > k = dim(Null(,B)). Since dim(Null(, 8 )) > d + 1 according to Lemma 2, we have 
Ua > dim(Null(, 8 ))/(i > (d + l)/d > 1. □ 


A simple but important fact suggested by Corollary 1 is that any localizable network must have at least two anchors. 
Similar conclusions have already been obtained in the existing studies for networks in the two-dimensional space 
[7-10]. But Corollary 1 also suggests another important fact, which has not been observed in the literature, that 
more anchors are required to ensure the localizability when dim(Null(, 8 )) increases. The quantity dim(Null(,B)) can 
be viewed as a measure of the “degree of bearing rigidity” because dim(Null(, 8 )) reaches the smallest value d -I- 1 
when the network is infinitesimally bearing rigid as shown in Lemma 2. As a result, the intuition behind the second 
fact is that more anchors are required to ensure the localizability when the network is “less” bearing rigid (i.e., 
dim(Null(, 8 )) is large). 

We next present another three localizability conditions, two of which are sufficient and the other is both necessary 
and sufficient. These conditions are important because they indicate the explicit connection between the localizability 
and infinitesimal bearing rigidity. Before presenting the conditions, we need to first define the notion of augmented 
network. 

Definition 2 j]Augmented Network) Given a network Gip) with Q = {V,S), denote by G{p) an augmented 
network with G — {V,S) where S = £ U {{i,j) ■ i,j G Va}- 

The augmented network G{p) is obtained from G{p) by connecting every pair of anchors. If the anchors are already 
connected in G (p) , then G(p) is the same as G (p) ■ It should be note^ that adding or deleting the edge between any pair 
of anchors only changes Baa but not Bjf. As a result, G{p) and G(p) have exactly the same Bff and hence they are 
localizable or nonlocalizable simultaneously. The next two sufficient conditions connect the notions of localizability 
and infinitesimal bearing rigidity. 

Corollary 2 When Ua > 2, ifG{p) is infinitesimally bearing rigid, then G{p) is localizable. 


Proof. We will first use TJreorem 2 to prove the localizability oi_G{p). Then the localizability oi_G{p) immediately 
follows because Gip) and Gip) have the saine localizability. Let B be the bearing Laplacian for Gip). Since Gip) is 
infinitesimally bearing rigid, we have Null(, 8 ) = span{l 0 IcItP} by Lemma 2. As a result, any infinitesimal bearing 
motion Sp G Null(,B) can be expressed as a linear combination of 1 0 and p. Since no two anchors collocate, there 
does not exist a linear combination of 1 0 and p leading to Spa = 0 if no > 2 . Then Gip) is localizable according 
to Theorem 2. □ 



Corollary 3 When Ua > 2, if Q{p) is infinitesimally hearing rigid, then Q{p) is localizable. 


Proof. Similar to Corollary 2. □ 


The intuition behind Corollary 3 is as follows. If a network is infinitesimally bearing rigid, then it can be uniquely 
determined up to a translation and a scaling factor by the bearings. Since the translational and scaling ambiguity can 
be further eliminated by the anchor constraints, the entire network can be fully determined and hence localizable. 
It is notable that Corollary 3 is more restrictive than Corollary 2 because it requires G{p) to be infinitesimally 
bearing rigid whereas Corollary 2 merely requires Q{p) to be. To illustrate, for each of the networks as shown in 
Figure 3(c)-(f), the augmented network Q{p) is infinitesimally bearing rigid but G{p) is not. Then, these networks 
can be concluded as localizable by Corollary 2. Finally, Corollary 2 can be viewed as a generalization of the result 
[10, Cor 10] which is applicable only to two-dimensional cases. 

As suggested by Corollary 2, the condition of the infinitesimal bearing rigidity of Gip) is sufficient to ensure the 
localizability of G (p) ■ An important yet unexplored problem is whether or not the condition is also necessary. In the 
case of Ua > 3, the condition is sufficient but not necessary. For example, for the network in Figure 3(g), Gip) is 
localizable but G (p) is not infinitesimally bearing rigid since the three anchors are collinear. However, in the case of 
Ua = 2, the condition is both necessary and sufficient as shown below. 

Theorem 3 When Ua = 2, a network G(jp) is localizable if and only if the augmented network G{p) is infinitesimal 
bearing rigid. 


Proof. The sufficiency_has already been proved in Corollary 2. We next prove the necessity by contradiction. Assume 
G{p) is localizable but Gip) is not infinitesimal bearing rigid. Then G(p) has a nontrivial infinitesimal bearing motion 
Sp which is not in span {1 0 Id,p}. Write 6p = [ 5 pJ, (*)]"'", where 6pi,6p2 € corresponds to the two anchors. 

Because the infinitesimal motion Sp preserves all the bearings including the bearing between pi and p2, we know that 
the vector Spi — 6p2 is parallel to Pi—p2- Asa result, there exists a nonzero scalar k such that Spi — 6p2 = fc(pi — ^ 2 )- 
Construct 


Sp = Sp + In ^ {kp2 — Sp2) — kp 


Spi 


kp2 - Sp2 


kpi 


0 

Sp2 

+ 

kp2 - Sp2 

- 

kp2 

= 

0 

A) _ 


(*) 


_(*)_ 


_(*)_ 


Since the first two entries of Sp' are zero, we know Sp' is an infinitesimal motion that only involves the followers. 
Thus, the network is not localizable by Theorem 2, which is a contradiction. □ 


5 Distributed Network Localization Protocols 

In this section, we propose and analyze a linear distributed protocol for bearing-based network localization in 
arbitrary dimensions. 

The global minimizer of the unconstrained optimization problem (5) can be obtained by the gradient decent protocol 

Pfi^) = J(p/) = -BffPfit) - BfaPa, (6) 

whose elementwise expression is 

Mt) = -Yl Pg^AMt) - Pjit)), * € V/. (7) 
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Fig. 4. The geometric interpretation of protocol (7). 

where Pg.. = Id — 9ij9ij- Note the neighbor of the follower i can be either a follower or an anchor. 

Several remarks for protocol (7) are given below. First, the protocol is distributed because the localization of Pi 
only requires {9ij}jeMi and In practical implementation, the bearings {9ij}jeJ9'i can be measured by 

a bearing-only sensor such as a camera and the estimates {Pj}jeN'i can be transmitted from the neighbors via 
wireless communication. All the bearings must be measured in a global reference frame. Second, the protocol has a 
clear geometric interpretation as shown in Figure 4. The term —Pg..{pi(t) — pj(t)) is the orthogonal projection of 
(pjit) — pi{t)) onto the orthogonal compliment of pij, and hence it acts to steer the estimate Pi{t) to align with the 
bearing measurement pij. Third, protocol (7) can be viewed as an extension of the protocol proposed in [10], which is 
applicable to networks in the two-dimensional space. Finally, those who are familiar with consensus problems might 
have noticed that protocol (6) has a similar expression as the well-known consensus protocol [18]. The difference is 
that in the consensus protocol, the weight for each edge is a positive scalar whereas in the localization protocol the 
weight for each edge is a positive semi-definite orthogonal projection matrix. 

The convergence of the protocol is characterized as below. 

Theorem 4 The distributed protocol (7) can globally localize the network Q (p) if and only if the network is localizable. 

Proof. When Bff is nonsingular (i.e., the network is localizable), the matrix —Bff is Hurwitz. As a result, the 
linear time-invariant system (6) is stable and the state converges to the steady state value —Bj^BfaPa which equals 
to the real follower location pf according to Lemma 3. When Bff is singular (i.e., the network is not localizable), 
the final estimate would depend on the initial estimate of the network location. □ 

5.1 Sensitivity Analysis 

Since the bearing measurements may be corrupted by errors in practice, it is meaningful to study the impact of 
constant measurement errors on the localization protocol (7). Denote the unit vector G as the measurement 
of gij. In the presence of bearing measurement errors, the localization protocol (6) becomes 

Pf(i) = - BfaPa, (8) 

where Bff and Bfa are obtained from Bff and Bfa by replacing gij with gij, respectively. The matrix Bff may not 
be symmetric since gij —gji in general. 

We next analyze two problems regarding (8). The first is when Bff is positive stable (i.e., all its eigenvalues have 
positive real parts) such that (8) is globally stable. If Bff is positive stable, the final estimate given by (8) is 

P*f = -BjjBfaPa. (9) 

The second problem is how large the localization error jjpj ~ P/|! is. To solve the two problems, define 

ABff^Bff-Bff, ABfa^Bf^-Bfa, 

as the perturbations of Bff and Bfa caused by the bearing measurement errors. Let 0ij G [0,7r] be the angle 
between gij and gij; that is gfjgtj = cos Oij. The angle Oij represents the inconsistency between gij and gij. This 
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representation is valid for arbitrary dimensions. Note 9ij ^ 0ji in general. Define the total bearing measurement 
error for the followers as 


ieVf j&Mi 

We next give lemmas to characterize the relationship between e and ^Bff, ^Bfa- 

Lemma 5 Denote by 6 G [0,7r] the angle between any two nonzero vectors x,y G (i.e., x^y = ||a;||||y|| cos6*J. 
Then \\Px — Py\\ = sin6l. 


Proof. See Appendix B. □ 

Lemma 6 For a network Q{p) with arbitrary hearing measurements {gij}{ij)^£, it always holds that |jA,Byj:|| < e 
and |jA,Bya|| A e/2. 


Proof. Denote ^Pgij — Pgij —Pgij,^ihj) € £■ K then follows from Lemma 5 that ||APg^^ || = sin^ij. Note [ABff]ii = 
^Pgij fo'' * ^ ''2/i [^Bff]ij = —APg.. for i GVf and j G AfiCi V/; and [ABff]ij = 0 otherwise. Then we have 

AP„.J| + llAP„.Ji < 




I1AP//II < J2i^Vf + Si6V/ 

2 II “ 2 sin0, 

J2ieVf ll^-^so II “ J2ieVf sin% = e/2. □ 


ij = e. Similarly, we have \\ABfa 


9ij I 


V/ 




A Z^ieVf SjeTVinVa ll^-^9o II - 


We now give a upper bound for the total bearing error e to ensure the positive stability of Bff. 

Theorem 5 Given a localizable network with Bff nonsingular, the matrix Bff is positive stable if the total hearing 
error e satisfies 


e < Aniin(^//)) 

where Ami„(P/y) is the minimum eigenvalue of Bff. 


( 10 ) 


Proof. Since IIAP//I1 < ebyLemma6,if (10) holds, we have ||AP//|| < Amin(^//) = l/lIPj^H, which further implies 
\\BjjABff\\ < IjP^^IIIIAP//II < 1. Thus the spectral radius p{Bj^ABff) < 1 and hence the matrix (/ + Bj^ABff) 
is nonsingular. As a result, Bff = Bff + ABff = Bff{I + Bj^ABff) is nonsingular. Since Bff is obtained by 
perturbing Bff and Bff is positive stable, the nonsingularity of Bff implies the positive stability. □ 

Theorem 5 suggests that a large Amin('8//) would give the network a large tolerance to bearing measurement errors. 

We now study the localization error ||pj ~P/||- An intuitive conclusion that can be immediately drawn from (9) and 
matrix perturbation theory is that the localization error would be sufficiently small when the bearing measurement 
errors are sufficiently small. We next give a specific upper bound on the localization error. 

Theorem 6 The estimate p*f = —BjjBfaPa given in (9) satisfies \\p*f — p/|| < ^ (ilboll + Ib/ll)- 


Proof. See Appendix C. □ 


In the last, we briefly discuss the impact of measurement errors in the anchors’ locations. Suppose the bearing 
measurements are accurate in this case. Then the final estimate given by protocol (7) becomes = —Bj^Bfa{pa + 
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yimj 


yimj 





Fig. 5. Simulation examples for the localization protocol (7). The bearing measurements are accurate for the example in 
(b)-(c), and inaccurate for the one in (d)-(e). The blue squares represent the anchors. The blue hollow dots and the green 
solid dots represent the true and estimated locations of the followers, respectively. 

/S.pa), where G denotes the anchor location error. Then the localization error is given by /S.pf = pj — pf = 
—BjjBfa^Pa, which indicates that the anchor location errors prorogate to the final localization error via a linear 
transformation. It is straightforward to show that a translational or scaling error in the anchor measurements would 
cause the same translational or scaling error in the localization of followers. 


5.2 Simulation Examples 


Two simulation examples are shown in Figure 5 to demonstrate the localization protocol (7). The network to be 
localized is a three-dimensional cubic network, which contains eight nodes and two of them are anchors and the other 
six are followers. The initial estimate, which is randomly generated, is given in Figure 5(a). For the first example in 
Figure 5(b)-(c), the bearing measurements are accurate and it can be seen that the estimate of the network location 
converges to the true value. For the second example in Figure 5(d)-(e), the bearing measurements are inaccurate. 
Specifically, the total bearing error is e = 2.77 and the final localization error equals 7.25 m. By comparing the 
two examples, it can be seen that when the bearings have measurement errors, the finally localized network would 
have localization errors. However, the final localized network can still be sufhciently close to the true network if the 
bearing errors are sufficiently small. In addition, for the second example, we have Amin = 0.59 < e. Although the 
condition in Theorem 5 is not satisfied, the matrix Bff is still positive stable which indicates that the condition in 
Theorem 5 may be conservative. 


6 Conclusions 

This paper studied the problem of bearing-based network localization in arbitrary dimensions. The first main con¬ 
tribution of this paper is to propose a variety of necessary and/or sufficient conditions for network localizability. The 
second main contribution is to propose and analyze a linear localization protocol. The results presented in this paper 
not only can be applied to solve the problem of sensor network localization but also provide a theoretical foundation 
for bearing-based multi-agent formation control [17,19-21]. 

In this paper, we assumed that the underlying graph is undirected. As we have explained, the localizability analysis 
is independent to whether or not the sensing graph is undirected because any directed graph can be converted to an 
undirected one without affecting the localizability analysis. But the convergence analysis of the proposed localization 
protocol relies on the assumption of undirected graphs. As observed in [17], a new notion termed bearing persistence 
emerges and makes the problem more complicated to analyze in the directed case. Distributed localization with 
directed interaction topologies is therefore a direction for future work. 
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A Preliminaries to Bearing Rigidity Theory 

For a network G{p), consider an oriented graph and express the edge vector and the bearing for the fcth directed edge in the 
oriented graph, respectively, as and = efc/||efe|| for k £ {l,...,m}. Define the bearing function Fb : R'*” —>■ R‘*"‘ as 
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Fb{p) — [gl, ■ ■ ■ , The bearing rigidity matrix is defined as the Jacobian of the bearing function, Rb{p) — dFB{p)/dp € 

Ijtimxdn important properties of the bearing rigidity matrix are given as below. 

Lemma 7 ([15]) For any network Q(p), the bearing rigidity matrix satisfies Rb = diag {Pg^/Wek\\) H, Rank(i?s) < dn — d—1 
and span{l ® Id,p} ^ Null(7?s). 

Let Sp be a variation of p. If Rb{p)Sp — 0, then 5p is called an infinitesimal hearing motion of Q(p). A network always has 
two kinds of trivial infinitesimal bearing motions: translation and scaling of the entire network. 

Definition 3 (Infinitesimal Bearing Rigidity) A network is infinitesimally bearing rigid if all the infinitesimal bearing 
motions are trivial. 

The necessary and sufficient conditions for infinitesimal bearing rigidity are summarized as below. 

Theorem 7 ([15]) For any network Q(p), the following statements are equivalent: 

(a) Q{p) is infinitesimally bearing rigid; 

(b) G{p) can be uniquely determined up to a translation and a scaling factor by the inter-neighbor bearings; 

(c) Rank(i?B) = dn — d — 1; 

(d) Null(i?s) = span{l (g) 

B Proof of Lemma 5 


Proof. Here we only prove the case of d = 3. Without loss of generality, assume x and y are two unit vectors satisfying 
ll®ll = lidll = 1- Then, we have Px = Id — , Py = Id — yy"^, and hence \\Px — Py\\ = Wxx"^ — yy'^W- There always exists 

an orthogonal matrix U £ such that the two vectors x and y can be orthogonally transformed to Ux = [1,0, ,0]"’’ and 

Uy = [cos d, sin d, 0]"^. Since the spectral norm is invariant to orthogonal matrices, we have 


P — P 

-I X -I y 


= \\U{xx'^ - yy^)U'^ 



1 0 



0 0 


= sind <5 

1, 


cos^ 9 
sin 6 cos 9 


cos 9 sin 9 
sirfi 9 


where Q = 


sind — cos 9 


— cos 9 — sin 9 
sin0||7|| = sind. □ 


It is easy to see Q^Q 


I2 and hence Q is an orthogonal matrix. Then, \\Px — Py 


sin0|]Q|| = 


C Proof of Theorem 6 


dl pf = —BjjBfaPa- Rewrite p} as p*f = —{Bff + ABff) ^{Bfa + ABfa)Pa. By [22, Eq. (25)], we have 
{Bff-\- ABff)~^ — Bfj — BJ^ ABff{I Bjj ABf f)~^ Bjj , substituting which into p*f gives p*f = —BjjBfaPa — BjjABfaPa + 
B7}ABff{I + Bj}ABff)-^Bj}ABfaPa + BjjABffil + Bj}ABff)-^Bj}BfaPa = Pf - {I + Bj}ABff)-^B7}ABfaPa + 




BflABff{I + BflABff) ^pf. It follows that 


Up} - p/ll < \\{I + Bj}ABff)-^Bj)ABfaPa\\ 

+ \\Bj}ABff{I + Bj}ABff)-^Pf\\ 

< II (/ + BjlABff)-^ II \\BJI II WABfx II ||p„ II 
+ ||I?7/||||AI?//||||(/ + S7/AI?//)-i||||p/|| 

= ||(/ + S7/AS,/)-^||||I?7/||(|lAI?/J|||p,,|| 

+ I|AS//|1||P/||) 
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Substituting ||AZ?//|| < e and |jA^?/a|| < e/2 as shown in Lemma 6, and |j(/ + BjjABff) ^|| < 1/(1— || |jAZ?//||) by 

[23, Lemma 2.3.3] into the above inequality gives 


\\P}-Pf\\ < 


mj}\mPa\\ + \\Pf\\)e 

l-\\Bjl\\e 


Substituting = 1/Amin(^//) completes the proof. □ 
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